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We reexamine the estimate of the neutron Electric Dipole Moment (NEDM) from chiral and QCD spectral sum rules (QSSR) 
approaches. In the former, we evaluate the pion mass corrections which are about 5% of the leading Log. results. However, 
the chiral estimate can be affected by the unknown value of the renormalizaton scale v. For QSSR, we analyze the effect 
of the nucleon interpolating currents on the existing predictions. We conclude that previous QSSR results are not obtained 
within the optimal choice of these operators, which lead to an overestimate of these results by about a factor 4. The weakest 
upper bound \Q\ < 2 x 10 -9 for the strong CP-violating angle is obtained from QSSR, while the strongest upper bound 
\6\ < 1.3 x 10 -10 comes from the chiral approach evaluated at the scale v = Mm- We also re-estimate the proton magnetic 
susceptibility, which is an important input in the QSSR estimate of the NEDM. 



1. Introduction 

The Lagrangian of Yang-Mills theory contains, in 
addition to the usual term, also a topological term: 



C = ~F^F^ -6q(x) 



(1) 



where q(x) is the topological charge density given by: 



q{x) 



32tt 2 



p»v = l e w Fpa . ( 2 ) 



The additional term violates the invariance under CP. 
This is called strong CP violation [I] to distinguish it 
from the CP violation present in the weak and elec- 
tromagnetic sector of the Standard Model. Experi- 
ments, however, do not show any violation of strong 
CP and require a very small value for \8\ < 2 x 10 -9 
as we shall see later on. In this paper we shall deter- 
mine the dependence of physical quantities on and 
study the processes that violate strong CP. For this 
purpose, we reevaluate the neutron electric dipole mo- 
ment (NEDM) which depends on the ttNN coupling 
which violates CP. 

2. Improved chiral estimate of the NEMD 

An elegant way of doing this is to use the low en- 
ergy effective Lagrangian of QCD that contains the 
fields of the pseudoscalar mesons and baryons instead 
of the original quarks and gluons. This is due to 
the fact that in the effective Lagrangian the effect 
of the axial U(l) anomaly is explicitly displayed and 
because of this the amplitudes for the hadronic pro- 
cesses can be easily computed. This Lagrangian can- 
not be explicitly derived from the fundamental QCD 
Lagrangian as in the CP^ -1 model (see e.g. Ref. |2I3) 
and References therein), but can only be constructed 
requiring that it has the same anomalous and non- 
anomalous symmetries of the fundamental QCD La- 
grangian. The expressions of the QCD effective La- 
grangian describing pseudoscalar mesons and includ- 



ing the U{\) anomaly are given in 
and the review in Ref. |2l3j . 
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• Estimate of the irNN couplings 
The ttNN couplings are defined asd: 

CkNN = VZttNt [i^g^NN + S-kMn) N , (3) 

where r are isospin Pauli matrices. 

The CP-conserving coupling is well measured : 

girNN = 13.4. (4) 

The CP- violating coupling g^NN can be obtained 
using an effective Lagrangian approach, though this 
approach for estimating the coupling can be question- 
able. It reads [7151: 



g-nNN 



m u m d 



— (m E - TO S ) X 

U{m u + md){m s - to) 

3m(TO£ — toa) 



2(m s — m)rriN 



(•5) 



with: 2to = (m u + md). We shall use the ChPT mass 
ratio [TS] : 



2to, 



''3 



m d ) 



24.4 ± 1.5 



(6) 



which is confirmed by the QSSR estimates of the quark 
mass absolute values in units of MeV [12ll3j : 

m u (2) = (2.8±0.2) , m d (2) = (5.1 ± 0.4) , (7) 

and: 

to s (2) = (96.3 ± 17.5) . (8) 

Using — 92.46 MeV, one can deduce: 

g„ NN = (0.0282 ± 0.0071)61 , (9) 

which is relatively small compared to g^NN- One can 
also note that the chiral correction is about 0.5% of 
the leading order expression given in Ref. [7J. We 
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2 We follow the normalizations of [7] but adding an overall \/2 
factor for charged pion fields. 
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have added an error of about 25% as a guess of the 
systematics of the method based on its deviation for 
predicting the value of the measured CP-conserving 
coupling Q 

• The NEDM from chiral and ChPT approaches 

To first order in 6, the neutron electric dipole mo- 
ment (NDEM) D N is given by: 

= T(n{ Pf )\J^0) i J d 4 xSC C p(x)\n( Pl )) 

= - {iD N ) u{ Pf )<T^k» l5 u( Pl ) + 0{k 2 ) , (10) 

where : k = p t — Pi is the photon momentum and: 

8C C p{x) = q(i-j b )Aq , (11) 

with A a 3-dimension hermitian matrix acting on the 
flavour space (u, d, s) and = (1/2) (7^ - 7^) . 
In order to extract Djy, we use the Gordon decompo- 
sition for the axial current: 



choice of the cut-off M^. One can also note the non- 
analytic terms y/a and a 3 / 2 . In the same way, Fig [5] 
gives: 

= i + ii±^/ arctan . ( 2 -«) 



arctan ■ 



V^4 
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a \ a 

— a)a > £ 



V( 4 - a ) a 



l-f^+^(3-loga) + T |™ 3 / 2 . (18) 



q ? + * 



p + fc 



r p + At 



u(p+k)j 5 a^k l/ u(p) 



-M(p+fc) 75 (2p+fc) MM (p).(12) 



The lowest order contribution to comes from the 
diagrams in Figs [1] and [2] Using the expression of the 
77777 vertex : 



Figure 1. 1st class of diagrams contributing to Dn- 



(tt-(p + fc)|J M |7r-(p)) =-(2p+ k) M + 0{k 2 ) 



(13) 



the evaluation of the previous diagrams can be ex- 
pressed as: 
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(14) 



where are integrals over Feynman parameters 

coming from Figs [T] and [5] Fig [T] gives: 



7* 1 



sc(l - y) 



1 ,1 

xdcc / dy 

" Jo x 2 (l - y) 2 + a 1 - x(l - y) 



log a 



(1 - a) + /(a) 



(15) 



with: 
/(«) 



a (3 — a) 

V a ( 4 - a ) 
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and where: 



a = {m^/M N y 



(16) 



(17) 



One can notice that the term logo appears naturally 
in the unexpanded full expression without an arbitrary 
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q ~> q + k ^7 + k "/T^ 9 ^ k 

k\< k 



Figure 2. 2nd class of diagrams contributing to Dn- 



• Results and discussions 

Adding the previous Feynman integrals 1^ and 1^ 
[ Eqs. dm) and ([18])] into the expression of NEDM in 
Eq. (|14p . the NEDM from chiral approach reads: 



D 



JV chiral 
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(19) 



The leading-log term agrees with the original result in 
Ref [7] . The cancellation of the constant terms have 
been also noticed in [16]. In addtion, we also have 
a cancellation of the a log a term, which implies small 
mass corrections (about 5%). Our result demonstrates 
the accuracy of the leading-log approximation used in 
Ref. 0. It also shows that the scale, appearing in the 
leading-log a term, is the mass of the nucleon but not 



3 



any arbitrary cut-off scale, because this term appears 
before the expansion in a [see Eq. Q15[)]. This leads to 
the prediction: 



AvUrai = (20 ± 5) x 1O" 3 GeV- 1 , 
= (40 ± 10) x 10" 17 6» cm . 



(20) 



However, within Chiral Perturbation Theory (ChPT), 
the addition of counterterms in the effective La- 
grangian induces a log(MAr/V) 2 -term and leads to the 
(renormalized) NEDM expression [16]: 



g-irNNg-rrNN 




Dn{v)\ ChPT 



where v is an arbitrary hadronic scale and k an un- 
known constant. This result indicates that, only the 
coefficient of the logm„- term is model independent. 

For a conservative estimate and for a model indepen- 
dent result, we keep only the logm T term and move, 
like in [H] , the scale v above m n from the value of the 
constituent quark mass, which we take to be about 
(Mat/3) (within a 30% error), to M N = 940 MeV. Us- 
ing the values of the parameters in Eqs. (TJ} and ©, 
one obtains, in units of the electric charge e: 



D N (v) = (2.7 - 23.8) x 10 3 6» GeV" 1 
= (5.4 - 47.6) x 10 17 6» cm , 



(22) 



where the range takes into account the assumed 25% 
systematic uncertainties for estimating g^NN and the 
assumed 30% one for the value of the quark constituent 
mass. The value in Eq. ([20)1 corresponds to the value 
v = Afjv and where the small chiral mass corrections 
have been included. One can notice that the width of 
the range depends strongly on the unknown value of v, 
and where a fine tuning is obtained for its low values. 

3. NEDM and Magnetic Susceptibility from 
QSSR 

Alternative estimates of the NEDM have been done 
using QCD spectral sum rules [T7]. We shall re- 
examine these results in this section and present some 
alternative new sum rules. 

The analysis is based on the baryon two-point corre- 
lator put in a background with nonzero 9 and electro- 
magnetic field F^ v : 

S(q 2 )\e, F=i [d 4 x e^(0|TA^)iV(0)|0)k F (23) 



(24) 



N is the nucleon operators which can be written in 
general as : 

N(x) = 2{(ipCj 5 ip) ip + biipCip)^} , 



where C is the charge conjugate; ip is the quark field; 
6 is (a priori) an arbitrary mixing between the two 



4 In order to help the reader, we use the same notations and 
normalizations as in Ref. | 17[ . which is the same as the ones 
in I19| . The normalization of the QCD correlator differs by a 
factor 8 from |20I21| . 



operators: b = in the non-relativistic limit, which 
is the choice used in different lattice calculations [IS] . 
Originally, Ioffe [TH] has used the choice b = — 1 in the 
1st QSSR applications to the nucleon channel, which 
he has justified for a better convergence of the QCD 
series m the OPE. In [20121] EL a more general analysis 
has been performed by letting b as a free parameter 
and then looks for the 6-value where the result is less 
sensitive to the variation of b. From the overall fit us- 
ing different form of the sum rules, the optimal result 
for the nucleon mass and residue have been obtained 
for b = —1/5, which can be qualitatively understood 
by taking the zero of the derivative in b when retaining 
only the lowest order contribution. 
For the analysis of the NEDM, Ref. [17] has used a 
value 6=1, which differs completely from the previ- 
ous choices. Within this choice, the authors impose 
the zero coefficient of an non-analytic mass- Log. in 
the next l/Q 2 -corrections (g 2 = — Q 2 ) to the lowest or- 
der contribution. However, the vanishing of the mass- 
Log, to leading order does not (a priori) guarantee 
the absence of this contribution to higher orders. In 
the following, we shall test the stability of the existing 
results versus the variation of b. 

• Expression of the two-point correlator 

For the forthcoming analysis, we shall keep the coef- 
ficient of the term {F^a^ , "f^q^} , in the Lorentz de- 
composition of the nucleon two-point correlator given 
in Eq. ([3]) (some alternative choices are also possible) . 
The QCD expression reads [T7] : 



s(Q 2 )\ 



2\|th 




C 2a ln^r-1 



21, 



xColn 



(25) 



v and v x correspond to an UV and a small mass arbi- 
trary subtractions; (ip^) is the quark condensate. The 
coefficient-functions are: 



C2a 



Co = {b+\f{Ae d -e u ) , 



4(6-l) 2 e d 1 + -(2k + 



C2b — — - 



(46 2 - 46 + 2)e d + (36 2 + 2b + l)e„ , (26) 



where e u , are the electric charge of the u and d 
quarks in units of e. %, £, £ are the magnetic suscep- 
tibilities of the QCD condensates which encode the 
electromagnetic field dependence of the two-point cor- 
relator. In units of e, they are defined as [2"2"j : 



(0|^ct m ^|0)| f 
3WyG^|0)| F 

g(o\h 5 ^G*m\F 



3 For reviews, see e.g. I13I18| . 



eqKF^iipip} , 
ieq^Fftv^i/i) , 



(27) 
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where g is the QCD coupling and G°^ v the gluon field 
strength. The size of these magnetic susceptibilities 
have been estimated in the literature using different 
methods. The values [23] : 

k~ -(0.34 ±0.10) , f- -(0.74 ±0.2) , (28) 

induce (a posteriori) small numerical corrections in 
the present analysis and will not be reconsidered. On 
the contrary, the dominant contribution comes from 
X, which is not known with a good accuracy: 

Nr. 



x[GeV~ 



-8 



Triangle anomaly [24] 



-8, -6 Laplace SR(LSR) 
-3.3 Light cone SR[25] , 



(29) 



which we shall reconsider later on. Our analysis gives 
the value in Eq. (T4"3"]) . which is in better agreement 
with the one obtained in [22124] , 

The phenomenological part of the correlator can be 
parametrized in the zero width approximation by: 



S(q 2 )\ 



exp 
9, F 



.4 



"QCD continuum" 



N 



(30) 



where An is the nucleon coupling to its corresponding 
current; A is an arbitrary coupling which parametrizes 
the single pole contributions; "QCD continuum" 
stands from the QCD smearing of excited state contri- 
butions and comes from the discontinuity of the QCD 
expression. 

• Estimate of the Magnetic susceptibility 

Considering the nucleon two-point correlator in Eq. 
© in presence of an external electromagnetic field, one 
can derive the following LSR (neglecting anomalous 
dimensions) for each invariants related to the struc- 
ture {a'uvP + po'nv) and iip^v — PuT/m) [H] (hereafter 
referred as IS): 



e u M\l 



Pi) 



3AF 



±e u (K~20-2e uX (M 2 -htf) 



1 r 

4 



~ X 2 -M 2 N /M 2 ( , a \ 

AnG Vm 2 p ) ' 



(31) 



and gluon condensates; M 2 = 0.8 GeV 2 [I9l20l2ll30l 
[T5] parametrizes the mixed quark-gluon condensate: 
(0|Vi ( T MI ,(A a /2)G^V|0) = M 2 (0|^|0); The QCD con- 
tinuum contribution from a threshold t c is quantified 

as: 



pn 



l+t c T 



{tcrf 



(*cT) 



(33) 



A pi Bp are the single proton pole coupling to the two- 
point correlator, while = 32tt 4 X 2 v is the coupling of 
the double proton pole; k, £ and x have been defined 
in Eq. (|27[) : /i p and /i° are the proton magnetic and 
anomalous magnetic moments. The sum rules for the 
neutron can be deduced from Eqs. ([3"T]) and ([32]) by 
the substitution: 

e d e u ; [i p , Hp > /i n ; A p , B p ► A n , B n .(34) 

Multiplying Eqs. (]31|) and ([32]) by e<j and each corre- 
sponding neutron sum rule by e u and then subtracting 
the proton and the corresponding neutron sum rules, 
IS deduce: 



P*p&d P>n^u 

4a 2 



Mtr 



B n e d )T 1 = 



Mi 



X 2 



'(el - el) . 



(35) 



In order to eliminate the single pole contribution, IS 
apply, to both sides of Eq. ([35]) . the operator: 



An 



d_ 



(36) 



Using the LSR expression of the proton coupling A at 
from the M/v-component F% of the two-point correla- 
tor [19120121] : 



\ 2 



2ar~ 



M N 

IS deduces for r _1 



(37) 



M 2 ■ 

m N . 



- (l - — 

3 I + 6 Ml 



4 / 2 a 
3 V + 3 M% 



2.96 



-1.93 



(38) 



in remarkable agreement with the experimental values: 



P 



aM 



N' 



2 Cd + 



(l~Po) + 



N< 



-M 2 N /M 2 



24M 4 



M 2 



B, 



(32) 



0.946 



where M 2 = 1/r is the SR variable; 
GeV is the proton mass; a = 47r 2 1 (0 1 -0-0 1 0} | ; b = 
47r(0|a s G 2 |0) ~ 0.87 GeV 4 [2~9"IT3] are the quark 



exp 
H'p 



2.79 , 



exp 



-1.91 , 



(39) 



despite the crude LO approximation used for getting 
these predictions. Including the OPE and anomalous 
dimension corrections, IS deduce, from Eq. (|35p . the 
predictions @: 



exp 



3.0 , 



exp 



-2.0 



(40) 



6 Some relations between the neutron anomalous magnetic and 
its electric dipole moments has been also derived using light- 
front QCD approach |31| . which will be interesting to check 
from some other methods. 
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and the correlated value: 



X ■ 



-8 GeV" 



(41) 



for £ ss k w . However, by examining the LSR used 
by the authors, we notice that these sum rules do not 
satisfy r-stability criteria such that (a priori) there is 
no good argument for extracting an optimal result. 
In order to check the previous result, we solve the two 
equations: 



X 



dr p 



and 



X 



dr 



CI 



(42) 



for each given value of t c . The functions Zy 7 have 
been defined in Eqs. ([3"T) and ([3"2"[) . We use the values 
of k and £ given in Eq. ([28]) but they do not affect 
much the conclusions like in the case of IS. The anal- 
ysis is shown in Fig. [3J where a common solution is 
reached at r = 0.4 GeV -2 , though the curves do not 
exhibit r-stability region. 




1 1.5 2 2.5 3 3.5 4 4.5 



Figure 3. Analysis of x using LSR for b=-l: a) r- 
dependence for t c = 3. GeV 2 ; b) independence of the 
common solution in (a) for t — A GeV -2 . 



Again like in the proton mass sum rule, the testability 
is reached around t c = 3 GeV 2 [21]. Taking as a con- 
servative estimate the range of values t c — 1.6 ~ 3. 
GeV 2 , where the lowest value corresponds to the be- 
ginning of r-stability for the determination of the pro- 
ton mass, we deduce the optimal estimate: 



X 



.5 ±1.0) GeV" 



(43) 



in good agreement with the IS previous value |22j and 
the one in |24j using the quark triangle anomaly. At 
t = 0.4 GeV" 2 where a common solution has been 
obtained, one expects a good convergence of the OPE 
and smaller effects of radiative corrections. We have 
used the choice b = — 1 which we expect to give a 
reliable result like in the previous cases of the proton 
mass and Dn discussed in the next section where the 
results are almost unchanged (within the errors) from 
b = — 1 to the optimal value -1/5 obtained in the case 
of the proton mass |20|21j . 



• Test of the LSR results of Ref. [17J for NEDM 

From the previous QCD and phenomenological ex- 
pressions of the two-point correlators, one can deduce 
the Borel/Laplace sum rule (LSR): 

C(t) = r(r)^J—(D N AT 



29m* 



32tp 



y N M N 



XC Q (1 - p ) 



- ln(r^ 2 ) +7b - 1 



a 



21, 



C 2 a X 



(44) 



where r 1 = M 2 is the LSR variable. Ref. [T7] uses 
either the value [19120121] : 

1 (1.05 ±0.1) GeV 6 , (45) 



\2 



(2n) 



or its LSR expression |19|2QI2l] from the q = 7 M g M 
part of the correlator: 

5 + 26 + 56 2 



(2tt) 4 A 



,-rM 



AT" 



04 



-T- 3 x 



(l-p 2 )+ir(a s G 2 )(l-p )T 2 



(46) 



However, due to its high-dependence on r, this sum 
rule is much affected by the form of the continuum 
such that we shall not consider it. Instead, we shall 
consider either the value in Eq. (|45[) . or the expression 
of the residue from the Mm part of the correlator: 

^2 



(2n) 4 \ N M N e 



-tMI 



(^)t- 2 (7 -2b -5b 2 ) 



(1 - Pl ) - 3A/ 2 r(l - 6 2 )(1 - p ) ,(47) 



\2 



which has a lesser dependence in r. 

We show the results in Fig[5]for the previous value of 
and using, as in Ref. 17j, (ipip) = -[0.225 GeV] 3 
for a better comparison. 

— For the choice b — 1 used in [T7] , one can see from 
Fig [3Ji that the optimal value is obtained at M 2 = 
t -1 ~ 0.5 GeV 2 , which is relatively low for justifying 
the convergence of both the OPE and the PT series 
in a s . Fig (4Jd shows, like in the case of the analysis 
of the proton mass, that the t c stability is reached at 
high-value of 3 GeV 2 but the estimate does not move 
much from the optimal value t c = 1.6 GeV 2 obtained 
in the proton mass sum rule [19120121] , In this case, 
one can deduce: 

0.36) GeV 4 , 



r 6=i 



- X (0.34. 



(48) 



which reproduces the result of |17j . Assuming, like 
in Ref. [T7], that the single pole contribution can be 
neglected (which we shall test in the next section), and 
using the value of \ = —5.7 GeV -2 used in [T7], one 
can deduce from Eqs. (f44[) and ([48)l : 

D N \ b=1 w 9 x 10 3 6* GeV -1 . (49) 

Though (almost) trivial, the previous test is necessary 
for calibrating our sum rule and for checking our inputs 
in the next analysis. 
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• Direct extraction of from a new LSR 

— By inspecting the LSR in Eq. (UK), one can also 
isolate Djy by working with the LSR: 



£ 2 (r) = —rC . 

We show the results of the analysis in Fig. O 



(57) 



Figure 4. Analysis of —r/x using LSR for b=l: a) 
M 2 -dependence for t c — 1.6 GeV 2 ; b) independence 
for M 2 = .5 GeV 2 . 



• New estimate of A and choice of the nucleon 
currents 

We shall reconsider the previous analysis by aban- 
doning the choice 6=1 for the nucleon current and by 
giving a new estimate of : 



ta{t) 



1 



\- N M N 



(50) 



This analysis is summarized in Fig. [5] where we have 
used the value of \m in Eq. (|45|) and the running con- 
densate value |12|13|26j : 

(0|^V|0)(M) ~ -[0.266 GeV] 3 ln 4/9 (M/A) , (51) 

with A ~ 350 MeV for 3-flavours. 

— One can notice that the result is optimal in 6 for: 



6~ -0.5 , 

and more conservatively in the range: 
- 1 < 6 < , 



(52) 



(53) 



which does not favour the choice 6 = 1 used in [17] . 
The previous range includes the conventional choices: 
-1 in [19] . -1/5 in [20 21 and the non-relativistic limit 
6 = used in lattice calculations 181. 



— The 2nd (important) assumption used in Ref. [T7] 
is the neglect of the contribution of the single pole 
controlled by the parameter A. By inspecting the LSR 
in Eq. ((44|) , one can isolate A by working with the new 
LSR: 



(54) 



One can see in Fig. [SJ: that for all ranges of 6, ta is 
much smaller than r justifying the assumption of [17] . 
At the optimal range of 6 values given previously, one 
can deduce : 



r/x = -0.09 , 



ta/x 



-0.02 . 



(55) 



Using the quark mass values in 12 13J and the previous 
value of x m Eq. (|43|) . one gets: 

D N ~ -0.11x0m* ~ (2.24±0.12)xl0" 3 (9 GeV" 1 .(56) 



One can notice that the sum rule stabilizes at r w 0.75 
GeV -2 , which is smaller than in the previous analysis, 
showing a better convergence of the OPE. 

— We also study the dependence of the result on 
the value of the IR scale v x . The optimal value corre- 
sponds to: 



(80 - 200) MeV 



(58) 



which has the size of a typical IR chiral scale (pion or 
constituent quark mass). 

Therefore, we deduce the optimal value: 



D N ~ -(0.10 ± 0.03 ±0.03)x6'm* 
~ (2.06 ± 0.08) x 10" 3 6» GeV" 1 



(59) 



where the first (resp. second) error comes from the lo- 
calization of the extremum in r (resp.) of the i c values 
which we take to be t c = (1.6 ± 0.2) GeV 2 around the 
value obtained from the proton mass sum rule |2QI21j . 

• NEDM results and systematics from QSSR 

The two results from the LSR are in good agreement 
and lead to the final estimate: 



AvIqssr - (2-15 ± 0.10) x lQ- 3 6 GeV~ 



(60) 



which agrees with the range spanned by the chiral and 
ChPT estimate in Eq. ([20]) . 

In order to analyze the systematic errors in the ap- 
proach, we estimate using vertex sum rules the well 
measured coupling g^NN — 13.4. We use the symmet- 
ric configuration of the hadronic vertex in [27 from 
which, one obtains the LSR: 



g-,NN{XNM N )(Uml)T 3 e- M « T 
■(m u + m d )(ipip} (2 



167T 



(61) 



Using the expression of X% in Eq. (|46|). one can deduce 
the LO sum rule: 

16^/^ Z2 + 86 + 86 2 



9-itNN 



M N \5 + 2b + 5b 2 
while the one in Eq. (|47| gives: 



(m u + m d ) /2 + 86+ fi ^ 2 

g-rrNN — : ~ T 



7 - 26 - 56 2 



(62) 



(63) 



Using a double pole dominance and neglecting the 
QCD continuum, Ref. [17] fixes from Eq. §2$ the 
nucleon operator mixing to be 6 = 1/t = 0.307 for the 
1st sum rule to reproduce the experimental value of 
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Figure 5. LSR analysis of —r/x (red : dash-dotted curve) and —ta/x (blue : continuous curve) : a) r = 1/M 2 - 
dependence for b — —0.5, t c — 1.6 GeV 2 ; b) t c dependence for b = —0.5, r = 2 GeV~ 2 ; c) 6-dependence for t c = 1.6 
GeV 2 , t = 2 GeV" 2 




100 150 200 250 



Figure 6. Analysis of using LSR a) r = 1/M 2 in GeV 2 -dependence for b = —0.5, t c = 1.6 GeV 2 ; b) b- 
dependence for r = 0.75 GcV~ 2 , t c = 1.6 GeV 2 ; c) v x -dependence of the optimal value in units of MeV. 



9-kNN- One can notice that this sum rule is not accu- 
rate due to its r-dependence. For r ~ 1 GeV -2 , the 
previous value of the quark mass evaluated at 1 GeV 
is (m u + m d )(l) ~ 10.9 MeV. Including the QCD con- 
tinuum contribution with t c — 1.6 GeV 2 , the 2nd sum 
rule gives: 



jNN 



8.9 



(64) 



We consider its deviation by 33% from the data as the 
systematic error of QSSR for this estimate 0- There- 
fore, we consider as a conservative estimate of Dn from 
QSSR: 



D N 



| QSSR ~ (2.15 ± 0.71) x 10- 3 6 GeV~ 
~ (4.24 ± 1.40) x 1Q- U e cm . 



(65) 



4. Constituent quark results 



For a qualitative comparison of the results from the 
chiral and QSSR approaches, we use a simple model 
where the constituent quark interacts with electromag- 
netic field. Then, one can write 221: 



(0|^ M ^|0)|j 



3 X^(0|^|0) 



d 4 pTr{S(p,M q )a^} , (66) 



'Note that a more precise estimate of g^NN including the con- 
tribution of the two first lowest quark and gluon condensate 
contributions is claimed in Ref. 28 from the 1st sum rule using 
a different configuration of the hadronic vertex. 



where S(p, M q ) is the quark propagator in presence of 
an electromagnetic field: 



S(p,M q 



i 



P - Mg 



-(0|^V|0)| 



1 



1 



P-M~^~p-M q '"p-M q 
^11 



2M q p-M q p — M q 



(67) 



where M q the quark constituent mass and fi q ~ 2 is 
its anomalous magnetic moment. Then, one can derive 
the relation l22l: 



x<o|W|o> = ^ml^lli + £ 



(68) 



Using this relation into the LSR expression of Dn, one 
can deduce in units of e: 



D 



N 



3M In 



M 



r- 2 e M » T 



1 



2 



0771* 



4.4 x 10~ 3 6» GeV- 1 w 8.7 x 10 17 6» cm, (69) 



where we have taken b = in the non-relativistic limit, 
t' 1 w M% w v 2 , and we have used M q w (200 - 300) 
MeV w v x . We assume that this crude non-relativistic 
approximation is known with an accuracy of about 
50%, which gives the final estimate: 

AT | Const quark « (8.7 ± 4.4) x 1O" 17 cm . (70) 
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This value can be compared with the one from more 
involved LSR analysis. This approximate formula may 
indicate that Dn is dominated by the non- analytic 
Log. contribution like in the case of the chiral esti- 
mate of [7] rederived in the previous section, but at 
the quark constituent level. 



5. Final range of the NEDM values 

The previous results from chiral and ChPT ap- 
proaches in Eqs. (HOI), from QSSR in Eq. JB5]) and 
from a naive quark constituent model in Eq. (|70[) 
are comparable. However, a more definite compari- 
son with the chiral and ChPT estimate requires a bet- 
ter control of the value of the renormalization scale v 
and an improved estimate of the CP violating ttNN 
coupling. Also, search for some other contributions 
beyond the standard OPE of QSSR like e.g. the one 
of the D = 2 dimension operator discussed [32] , may 
be required. 

Combining these previous results with the present 
experimental upper limit (in units of e) [33] : 



^|cx P <6.3x l(T 2b cm, 

one can deduce in units of 10 -10 : 



(71) 



9 < (1.6 ±0.4) [Chiral] : v = M N 

< (1.3 ~ 11.7) [ChPT] : M N /3 < v < M N 

< (6.9 ± 3.5) [Constituent quark] 

< (14.9 ±4.9) [QSSR] . (72) 

These results indicate that the weakest upper bound 
comes from QSSR, while the strongest upper bound 
comes from the chiral estimate evaluated at the scale 
v = Mn- Present lattice calculations are at an early 
stage 34 and may narrow the previous range of values 
in the future. 
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